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^ ' Abstract 

en ■ Modular cosmology is plagued by overproduction of unwanted relics, gravitinos and 

^~^ ', especially moduli, at relatively low energy scales. Thermal inflation provides a com- 

f— N I pelling solution to this moduli problem, but invalidates most baryogenesis scenarios. 

■rj" ' We propose a simple model in which the MSSM plus neutrino mass term {LHy)'^ is 
supplemented by a minimal flaton sector to drive the thermal inflation, and make two 



i-S^ \ crucial assumptions: the flaton vacuum expectation value generates the /i-term of the 

II ■ MSSM and ?n,|^ + m|^ < 0. The second assumption is particularly interesting in that it 

Oh. violates a well known constraint, implying that there exists a nearby deep non-MSSM 

(— I \ vacuum, and provides a clear signature of our model which can be tested at future 

particle accelerators. We show that our model leads to thermal inflation followed by 



Affleck-Dine leptogenensis along the LH^ flat direction. A key feature of our leptoge- 
nesis scenario is that the HuH^ flat direction is also induced to temporarily acquire a 
C^ ' large value, playing a crucial role in the leptogenesis, as well as dynamically shielding 

the field configuration from the deep non-MSSM minimum, ensuring that the fields 
relax into our MSSM vacuum. 



1 Introduction 

Among the many unresolved problems in cosmology, the origin of the matter-antimatter 
asymmetry, or baryon asymmetry, is the one we can most easily recognize because we could 
not exist without it. According to recent data from the WMAP experiment P], the asym- 
metry is 

— = (6.5 ±0.4) X 10"^°. (1) 

Since the time, about forty years ago, when Sakharov showed that three conditions must be 
satisfied to generate this asymmetry J2] , many authors have tried to explain this asymmetry 
with many different models. But we still don't have any idea which, if any, is the right 
one because baryogenesis usually depends on unknown physics beyond the reach of current 
experiments. 

Besides producing necessary relics such as the baryon asymmetry, a successful cosmolog- 
ical history must also avoid or dilute unwanted cosmological relics such as gravitinos and 
moduli 111 which can destroy the successful predictions of Big-Bang nucleosynthesis or 
over-close the universe, depending on the scale of supersymmetry breaking. In gravity medi- 
ated supersymmetry breaking scenarios, these particles have masses of order mEw ~ 10^ to 
10^ GeV and Plank scale suppressed couplings. This indicates their long life time and pos- 
sibility of disturbing Big Bang nucleosynthesis. Both particles may be thermally produced 
in dangerous abundances in the early universe, but such production might be avoided if the 
reheat temperature after inflation was sufficiently low. More seriously, moduli have another 
more dangerous source, production of coherent moduli oscillations with amplitude of order 
of the Plank scale when the Hubble parameter drops below the vacuum mass of the moduli 
|3|. In the case of gravity mediated supersymmetry breaking scenarios, this occurs at an 
energy scale of the order of the intermediate scale, A/mEw^Pi ~ 10^° to 10^^ GeV, i.e. after, 
or at best at, the end of the usual high energy inflation. Thus, we face a disaster unless 
there is a proper amount of late time entropy release. 

Fortunately, thermal inflation |ni|7| provides a compelling solution to this moduli problem. 
It is expected to occur after the usual high energy inflation, at an energy scale between the 
intermediate scale and the electroweak scale, releasing enough entropy to dilute the unwanted 
relics. It is driven by the potential energy of a 'flaton' field which is held at the origin by 
thermal effects. When the temperature drops sufficiently, the flaton rolls away from the 
origin, ending the thermal inflation. The flaton then oscillates about its vacuum expectation 
value of the order of 10^° to 10^^ GeV till its decay completes at a temperature of order GeV. 

There are several categories of baryogenesis models and it is important to consider them in 
the context of the above discussion of a complete and consistent cosmological history. One is 
baryogenesis by out of equilibrium decay of heavy particles. For example, GUT baryogenesis 
|H| uses particles with mass of the order of the GUT scale, Mqut ~ 10^^ GeV. Right-handed 
neutrino leptogenesis [9^ uses decay of heavy right-handed neutrinos of mass 10^^ to 10^^ 
GeV, and anomalous B + L violating electroweak processes mediated by sphalerons [TD| to 
convert the generated lepton asymmetry to baryon asymmetry. However, in both cases, very 
high energy scales are required to produce the heavy particles. This means the generated 
asymmetry would be diluted to negligible amounts by thermal inflation or any other form of 
entropy production required to dilute the unwanted relics that are produced at lower energy 



scales. 

Another category is electroweak baryogenesis [TT]. It uses the electroweak phase transi- 
tion and sphalerons. The attractive point of this mechanism is that it works at a low energy 
scale, the electroweak scale, so it is rather independent of unknown high energy physics. How- 
ever, it seems that the window in the Minimal Supersymmetric Standard Model (MSSM) for 
this mechanism is small J21- Moreover, assuming thermal inflation, the reheat temperature 
after flaton decay is too low even for this mechanism. 

The last major category, Affleck-Dine (AD) baryogenesis ^Hj is distinctive from the 
previous two categories due to its efficiency. Using coherent fields along MSSM fiat directions, 
it can generate huge asymmetries which can survive substantial dilution. Indeed, in the case 
of gauge mediated supersymmetry breaking, it can even produce a large enough asymmetry 
to survive the full entropy production of the thermal infiation needed to dilute the moduli 
to acceptable levels [T3|. However, especially in the case of gauge mediated supersymmetry 
breaking, the AD condensate can fragment into Q-balls fSlCEl; though it is not clear how 
serious a problem this is as the absorbed baryon asymmetry may be released through the 
decay of the Q-ball ^Tj. In the case of gravity mediated supersymmetry breaking scenarios, 
the required entropy production is greater, so that not even AD baryogenesis can survive 
the full entropy production of the thermal infiation. 

Thus, if we assume supersymmetry breaking is mediated by gravitational strength in- 
teractions and thermal infiation solves the moduli problem, all the baryogenesis models 
suggested so far do not work. Therefore, a new scenario of baryogenesis is required under 
these assumptions. This is our motivation. The new baryogenesis model must satisfy two 
requirements. The first is that, as it cannot work either before the thermal infiation or after 
the fiaton decay, it must work around the end of the thermal infiation. The other is that it 
must generate a large asymmetry which can survive the dilution by the entropy released in 
the fiaton decay. In this paper, we suggest such a model. It is similar to an earlier model by 
one of us [TH], but is based on a more minimal model and crucially relies on some dynamics 
neglected in [T5] . 

We begin our discussion by explaining the construction of our model in Section [21 In 
Section ini we describe our baryogenesis scenario. In Section |1[ we discuss the constraints on 
parameters and the stability of dangerous field directions. In Section we show the results 
of numerical simulations of the field dynamics. In Section we discuss some of the issues 
involved in determining whether our generated asymmetry is preserved. In Section [71 we 
conclude and comment about future work. 

2 The Model 

2.1 The super potential 

The MSSM superpotential is 

where i is the generation index, a is the SU(2) index and a is the SU(3) index. 



Neutrino oscillations |T3 120] imply nonzero neutrino masses. So we include the following 
Majorana neutrino mass term 

-XlhabLiH^ecdL^jH^. (3) 

We expect thermal inflation as a solution of the moduli problem. So we need a fiaton 
sector. For simplicity, we consider only a minimal fiaton sector^ consisting of a single field 
(j) with self-interaction 

^A^0l (4) 

This naturally stabilizes the fiaton field at a vacuum expectation value of the order of 10^'' 
to 10^^ GeV. The fiaton must also have some interactions with other fields, for example a 
coupling of the form A^0XX) fo provide the finite temperature potential which holds it at the 
origin during the thermal inflation, x ^i^nd x may carry MSSM charges but, as they become 
very heavy at the end of thermal inflation, are not MSSM fields. They should form complete 
SU(5) representations to preserve gauge coupling unification, but would still be expected 
to leave an imprint on the renormalization of MSSM parameters from the GUT scale. As 
these fields and couplings do not affect our dynamics beyond holding the fiaton at the origin 
during the thermal inflation, we do not specify them explicitly. 

The origin of /i in Eq. (j2)) is unknown. We assume it is absent initially and generated 
dynamically by the vacuum expectation value of the fiaton field via the interaction 

X,4>^ea,H:H', (5) 

as was suggested long ago in Ref. [21]. In our vacuum, this coupling gives 

which has the correct size for /z under the assumption 

lA^I ~ lA^I (7) 

and fiaton mass of electroweak scale. This coupling also gives a reheat temperature after 
fiaton decay high enough for thermalization and standard freeze out of neutralino dark 
matter jTS]. 

Given the structure of our model, we could also expect a term 

-\h [tabH^H^^) (8) 



2 



with 



However, we assume 



h\ ~ |A^| ~ lA^I . (9) 



|A.| > |A<^| (10) 



^One interesting alternative would be to embed a global U(l) Peccei-Quinn symmetry j221 in the fiaton 
sector. Our scenario would work similarly in that case, though there is a danger of overproduction of axions 
|2S! by the fiaton decay. 



to give correct values for the neutrino masses and an appropriate scale for thermal inflation. 
As we shall see, this means that \XH€abH!;^Hl^\ ^ w^ew throughout our dynamics, so such a 
term is negligible. 

Thus, our superpotential is 



+ X,<phabH:H', + \x^ea,LtHte,,^Ht + \x^ct^\ (11) 

We could use a Z4 symmetry with charges 

Z4{/f„, Hd, Q, u, d, L, e, 0} = {0, 2, 2, 2, 0, 0, 2, 1} (12) 

plus the usual i?-parity, to enforce the above form for the superpotential. This Z4 symmetry 
could be a gauge symmetry if it is anomaly free ^^. It can be made anomaly free by adding 
an extra two SU(5) 10 representations and two SU(5) 10 representations, with coupling 
A^^lOjlOj where i = 1,2 and j = 1,...,5. Some such coupling would in any case be 
necessary for the thermal inflation as already mentioned. Gauging the Z4 symmetry is one 
way to avoid the potential domain wall problem |2S] caused by the flaton potential's four 
minima. We assume we are free from that problem. 

2.2 Our ansatz 

In our vacuum, all the MSSM supersymmetric flat directions in fleld space are stable at 
the origin. However, during thermal inflation the flaton is at the origin, and so there is 
no /i-term. This may cause some supersymmetric flat directions involving Hu or Hd to be 
unstable during thermal inflation. These potentially unstable flat directions are HuH^ and 
LiHu- Thus, we assume that the only flelds which get nonzero expectation values are Hu-, 
Hii, Li and 0. We will check the consistency of this ansatz in Section 0] 

For simplicity, we will truncate to a single generation. This may be essentially correct 
if all the matrices of leptonic sector are simultaneously diagonalizable, while otherwise will 
hopefully still capture the essential points. 

Gauge flxing and imposing the D-term constraints for these flat directions, we can pa- 
rameterize them as 

//.-C;"). i/.=(,°J. L=(»). *=^ (13) 

with the remaining D-term constraint 

^ = |/ij'-|/id|'-|f = (14) 

and corresponding gauge degree of freedom. Integrating out the remaining gauge fleld and 
choosing the gauge A^ = Q gives 

Jm = ^ Kd.K - Kd.hl) - i {h*ad^hd - hddX) - ^ i^*d,l - ld,n = 0. (15) 



Thus, the superpotential reduces to 

W 
and the F and D-terni parts of the potential to 



Vp 



K<fKhd + \\J^hl + h -^^ 



2|2 



A/,0 hd + \yl hu\ + \\^,<^ K\ + \\^hu\ + V-^iJi^^uhd + A<^(/)' 



/,3|2 



|2\2 



where g^ = {gf + 5'|)/4. The vacuum supersymmetry breaking potential is 

Ksb 



(16) 

(17) 
(18) 
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Vo + mjjJhu\ +mH\hd\ +m 



l\l? + <W 



+ 



A^X^(P'Khd + -A^A^/^/i^ + -A, 



+ c.c. 



(19) 



where Vq is adjusted to give zero cosmological constant, and all the other supersymmetry 
breaking parameters are of the order of the electroweak scale. 

Again for simplicity, we neglect any field dependent renormalization of the supersymme- 
try breaking parameters. Note that taking into account the dependent renormalization of 
the supersymmetry breaking parameters will be important for accurate quantitative compar- 
ison of our constraints on the supersymmetry breaking parameters near the end of thermal 
inflation with those parameters in our vacuum. 

Thus, the whole potential for these fields is given by 



V = Ksb + Vf + Vd. 



(20) 



3 Baryogenesis 

We are now in a position to describe how baryogenesis can occur in our model. 



3.1 Before the end of thermal inflation 

At high temperature, the thermal masses hold all the fields at the origin. So, we take our 
initial condition as thermal inflation with all fields having zero expectation value. We assume 
m?, < and m\ + rn^u^ < so that the fiaton and the supersymmetric flat direction LHu 
become unstable as the temperature drops during thermal inflation. We assume LH^ rolls 
away first. Once a field rolls away, it decouples from the thermal bath, and we neglect any 
remaining thermal potential. We will discuss the wider implications of our rather non-trivial 
assumption that m\ + rn\j^ < in Section HJ 

While the fiaton is still held at the origin, the potential for / and hu is 



V = Vo + mi^JKr + mi\l\' + 

|2 



-A^XJ^hl + c.c. 

2\2 



+ \Xj'K\' + \XJhl\' + ^gmK\'-\l\') 



(21) 



which has miniinuiii at 

^.A,'2^2_ I, X ,2,2 

and 



A,XJ'hi = -\A,XJ'hi\ (22) 



\A,\ + J\A,\^-Q{ml + mlJ (^2 ^2 ) 



" ^ 6|A,| V ■ ^ ' 

For simphcity, we assume / and hu settle down to this minimum before or hd start 
rolhng away from the origin. Note that substituting Eqs. (|^ and (j^ into Eq. (fTHj) gives a 
contribution — (m| — m^^)/2 to /i^'s mass squared. We assume that this does not destabihze 
hd- A sufficient condition for this is 

"^k + "^L > 2 (""' + "'^") (2^) 

though it is probably not necessary. 

3.2 At the end of thermal inflation 

Thermal inflation ends when rolls away from the origin. We assumed |Ai,| ^ \\^\ to give 
correct values for the neutrino masses and an appropriate scale for thermal inflation. This 
condition also ensures that the flaton dynamics is dominant as the flaton gets a larger value 
than hui hd or /. Thus, the flaton can be regarded in isolation, neglecting any back reaction 
from the other fields. See Eq. ()56|) . The pure fiaton part of the potential is 



^Ti = K) + <i0r + 



-A^X^cj)^ + c.c. 



/,3|2 



A^0T (26) 



which has minima at = 00 plus Z4 symmetric points, where 



^0 
and 



^ - -\A^\^(t)l\ (27) 



-12m2 + |AJ2 +\A^ 
^0? 



6A, 



(2J 



From Eqs. (fTH|) . (fTTjl and (jIHI), once the fiaton settles down, the MSSM parameters /i and B 
will be given by 

/i = \<Pl (29) 



and 






The remainder of the potential, 



1 
2 



+ 



A^X^cjy'^Khd + -A^Kfhl + 2\l(f)*^X^(f)Khd + c.c. 



+ \\f,(f)'^hd + \j'^hu\ + |A^0^/i„| +\XJhl\ + \2X^(t)huhdf 

+ l9'{\K\'-\hd\'-\in' (31) 

determines the dynamics of hu, h^ and /. As (p rolls away from the origin, the A-term 
A^X^cfy^huhd + c.c. and the cross term Xll*'^h*^\^(f)^hd + c.c. shift the minimum of hd to 



'"a — 2 1 2 I 1 2 

"^H, - 2"^i + 2"^H„ 

The above equation is valid only for small 0. As and hd become larger, more terms become 
important and Eq. (jH^ only gives a rough estimate. This non-zero value of hd generated by 
is the key to generating a lepton asymmetry in our model, and is also crucial to maintain 
the stability of our ansatz as we will see in Section W?]\ This was missed in Ref. |18j . 

As approaches its minimum, the remaining terms become important. The A-term 
A^A<^0'^+c.c. determines the phase of 0"^, and the A-term A^X^cfy^huhd+c-c. and the cross term 
2X*^(j)*'^X^(j)huhd + c.c. determine the phase of (f^h^hd. The cross term A*0*^/i^A,^/^/i„ + c.c. 
rotates the phase of Ihu, which was initially determined by the A-term A^XJ'^hf^ + c.c, and 
at the same time \X^(j)'^hu\ gives Ih^ a positive mass squared bringing it back in towards 
the origin.^ The angular momentum of Ih^ generated in this way is our lepton number 
asymmetry. Note that the cross term X*(j)*'^h'^XJ^hu + c.c. must be larger than, or at least 
comparable to, A^XJ'^h'^ + c.c. to give an effective angular kick to the phase of Ihu. 

After this, the dynamics becomes somewhat complicated. We will analyze it numerically 
in Section El and discuss some of the relevant issues in Section IHl 

4 Constraints and Stability 

4.1 Summary of constraints on the parameters 

We assumed 

mj < (33) 

as is required for thermal inflation, and 

mi + mi^<0 (34) 

so that LHu becomes unstable near the end of the thermal inflation, allowing the possibility 
for leptogenesis. We will discuss the consistency and implications of Eq. (jH^ in Section 1^^ 

^In the case that the flaton sector contains the axion, terms equivalent to the A-term A^X^cj)'^ + c.c. and 
the cross term 2A!,(^* X^cjih^hd + c.c. may be absent. However, we see that the mechanism still works. 



We require any combination of LH^ and H^Ha to be stable at the origin. This gives the 
following constraints 25 . 

<„+<, + 2|/ir>2|i?/i| (35) 

if ^TT'L ^ ^TT'H + l/^P ~ l-^/^l SO that pure HuHd is the relevant direction, and 

«, + l/^P - ml) [ml + ml^ + l/ip) > \B^^\^ (36) 

if 'Til < m^j + l/ip — \B^\ so that some combination of LHu and HuHd is the relevant 
direction. If i? = 0, Eq. ()36|) reduces to the stability constraint along the pure LH^ direction 

< + "^k + l/^l'>0- (37) 

In our model, /i and B are given by Eqs. (|^ and (IHUj). Note that Eqs. (jHSj) and (|H^ are the 
well known MSSM constraints to avoid instability at large field values along a direction with 
nonzero L, Hu or H^ 25j. Our inclusion of the neutrino mass term does not make things 
more unstable. 

In addition to requiring that the origin is stable, we also have to worry whether the fields 
may get trapped in a minimum generated by the A-term A^XJ'^h'^ + c.c. It is not easy to 
determine explicit constraints that ensure that such a minimum does not exist for a general 
combination of LHu and HuHd directions, but the constraint for the pure LHu direction is 

< + '^k + l/^l'>l^-lV6. (38) 

For the general case, we can check whether the A-term minima exist or not for a given set of 
parameters. On the other hand, if the fields evade the minima dynamically, the constraints 
are not needed. But it is not clear in our model. 
Finally, we require 

(m|„ + l/ip) {ml^ + l/xp) < iB/xp (39) 

for correct electroweak symmetry breaking in our vacuum. 

4.2 Stability of our ansatz 

So far, we have analyzed our model and its dynamical behavior within the context of our 
ansatz that the only relevant fields are Hu, Hd, L and 0, assuming all other fields are zero. 
However, our key constraint, Eq. (|34j) . 

ml + mj,^<0 (40) 

violates the well known stability constraint for large field values along the directions of 
nonzero {Hu,L,Q,d) or {Hu,L,e) [211 123 • This is dangerous because it means that very 
deep non-MSSM minima exist in our model. Does this mean our assumption is invalid? 
The answer will be no if we live in our vacuum and its life time is longer than the age of 
our universe. Then two questions arise. The first is whether our field configuration will 
dynamically settle down into our vacuum, as opposed to the deeper non-MSSM one. If the 
answer to this question is yes, so that we live in a false vacuum, the second question is what 
is the life time of our vacuum. To answer the second question first, the MSSM has a large 



enough region of parameter space in which our vacuum is cosmologically stable so that we 
are safe fI7\ . 

Returning to the first question, one of the dangerous directions involves Qj" and d^. For 
simplicity, we restrict to a single generation. We can parameterize this single generation as 





Q 



V2 



d 



1 



V2 




The superpotential becomes 



w = - 










2/id + Xf,(f)^Khd + -XJ^hl + - 



(41) 



(42) 



and the potential 



V 



T/ I 2 II, |2 I 2 I I, |2 I -"^ / 2 I 2\ I |2 I 2|;|2 , 

Vo + mHjhu\ +mj^Jhd\ + - [mQ + m^) \q\ +m^|/| + 



m, 



+ 



1 11' 

-AqX^q^hd + A^X^(t)^Khd + -A^Xj'^hl + -A^X^cf)'^ + c.c. 



+ \X^(t)'^hd + Xj'^hu\ + 



Xdq^ + X^cty^K 



+ \2X^(t)Khd + X^(l)''\ +-g [\Kv + -\q 



+ \Xdhdq\ + \XJh^ 

2 



2|2 



ITi |2 |/|2 



(43) 



As rolls away from the origin at the end of thermal inflation, there is a danger that the term 
|Arfg^/2 + X^j,(f?hu\ may destabilize q rather than bring hu back in to the origin. However, as 
discussed in Section 1^^ the A-term AnXa4>^huhd + c.c. and cross term Xl,l*'^h*^Xn4P'hd + c.c. 
induce a non-zero value for hd, which tends to stabilize q via the term \Xdhdq\ . Thus, 
ignoring finite temperature effects which will also help to hold q at zero, q will be stable if 

|/p) >0 (44) 



{rriQ + ml) - |AA 



^ ^ a;0*2/^:/.s 



\hd\^ 



+ \Xdhdf' + -5' 



2 ''l/i |2 17, |2 



The D-term contribution is not expected to be large because the field configuration will 
follow a direction close to D-flat, and in any case would be expected to be positive, and the 
second term is not expected to dominate as Eq. ()32|) makes it roughly the geometric mean 
between the first and third terms. Thus one expects that Eq. (j44p will be satisfied, i.e. q is 
stable. We confirm this estimate using numerical simulations in the next Section. 
For the dangerous lepton direction 



V2 
I 



V 



1 



v^ 



we can apply the same argument as above with the substitutions 

q^v, Xd^ Ae, (mL m^) -^ {m\, ml). 



(45) 
(46) 
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5 Numerical Simulation of the Dynamics 

In this section we use numerical simulations of the homogeneous field dynamics to see the 
generation of the lepton asymmetry and check the stability of our ansatz. 

For our initial conditions, we take (p = h^ = 0, and / and h^ at the minimum of their 
potential. We take q slightly displaced from the origin to test its stability. We then displace 
(j) slightly from the origin at various angles, and follow the dynamics. To crudely mimic the 
transfer of energy from the initial homogeneous mode of to gradient energy of 0, which 
causes the homogeneous mode of to settle towards its minima, we add artificial damping 
acting on 0. Note that Hubble damping is negligible. 

There are various phases in our potential, but only certain combinations are physical. To 
make things more explicit, by field rotations we can choose 

arg (AxXx) = 0. (47) 

The remaining physical phases are arg (— A* A,y) and arg (— A* A^), with CP conserving values 
being and tt. 

We define the rescaled variables 

2 A 

i = mt, mj, = ^ = 0{l), Ax = — = 0{1), (48) 

" Mad ' ' Mad ' ^ Mad ' Mad ' '^ Mti ' ^ ' 

K = ^^^ = 0{1), A. = ^^^^ = 0(1) , A, = ^^ = 0(1), (50) 

rn m rn 

A^ = ^^^^^ = 0(10-1°) to 0(10-2) , (51) 

m 

gMAD ,^Mnb\ . A^/in7N ^ A^Mad 



C(IO^) to O(IO^) , \d = = C(l) to C(IO^) (52) 



m rn 



where 



M^ /\/f2 

^° - C(10-^) to 0(10-2) , -^ = C(l) to O(IO^). (53) 



m m 

The potential can be divided into two parts 



V = VTi{(p) + VAB{K,hd,q,l,(l)) (54) 

as in Eqs. pHjl and (jSH), and rescaled as 



T> _ ^TI f. _ VaD /p.p.N 



The rescaled Lagrangian is 



^ = zAfT (56) 



m^Mli 



i^TI^0 +^TI + 



Mad'' 



M 



TI 



^AD [dhu, dhd, dq, dVj + Vad (hu, hd, q, l, 0j 

(57) 
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where the K^s are the kinetic terms. Note that Vad only has a small effect on the dynamics 
of (j) because 

^] = O(lO-iO) to 0(10-2). (58) 

Mti / 

The equation of motion for (p is given by 

^ + r-^ + ^r^ = 59 

dt^ dt d(f)* 

with r our artificial damping. We do not include any damping for the AD fields as it is less 
clear what form it takes. 

Our numerical simulations show that q is stable. Thus the dynamics ensures that we do 
not fall into a deeper non-MSSM minimum, and our ansatz and key assumption, Eq. (|H^. 
are consistent. 

Some of our numerical results are shown in Figures d El and El 

6 Preserving the Lepton Asymmetry 

In the previous sections we showed that a lepton asymmetry can be generated in our model 
at the end of thermal inflation. However, there is a danger that this lepton number will be 
washed out by the continued action of the lepton violating operators. To avoid this, we need 
the amplitude of the lepton field to decrease, taking it into the conservative region of the 
potential. Therefore, some form of damping for the homogeneous mode of the lepton field 
is required to preserve the asymmetry. Note that Hubble damping, which is usually used in 
AD baryogenesis scenarios, is negligible. 

There are several sources of damping. The first we can expect is the transfer of energy 
from the homogeneous mode to inhomogeneous modes, i.e. the build up of gradient energy, 
which is often called preheating. It is clear that a substantial amount of energy would be 
transferred in this way in the first few oscillations, but it is not clear how long it would take 
this process to reduce the homogeneous field amplitude by a substantial factor. Friction 
induced by the thermal bath and decay are other sources of damping. We don't attempt to 
analyze these damping processes in this paper. 

Although we do not understand the details of the decay of the AD fields, we do know that 
it will complete while the flaton energy is still dominant, and will release enough thermal 
energy to restore the electroweak symmetry. Assuming our lepton asymmetry survives until 
this stage, it will be converted to a baryon asymmetry by the usual electroweak sphaleron 
processes. Finally, as the temperature drops to the GeV scale, the fiaton decay completes, 
releasing substantial entropy which dilutes the baryon asymmetry. 

Thus, the final baryon asymmetry is given by 

s n^ m^ Had n^ rn^ riAD \ Mti / rn^ ' 
As in Eq. (J3H|1 . we expect 

^^°^ ~ 10-10 to 10-2, (61) 



Mti 
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Figure 1: \(f)\ as a function of time, for m^, = 
arg0i = 37r/20, and artificial damping F = 0.2. 
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Figure 2: Dynamics of /, for the same parameters as Figure^ and rn|^ = —2.5, m^ = 2.4, 
ml = 1.0, \A^\ = 1.0, \A^\ = 1.0, |A^| = 1.0, |A^| = 2.5, g = 10^ arg (-A* A^) = tt + 7r/120, 
and arg (— A* A,^) = vr+vr/SOO. The correct value oi g is g ~ 10^, but the value we used is large 
enough to make little difference and is easier to handle numerically. Angular momentum 
is generated, but, as we did not include damping for the AD fields, we do not expect it to 
survive in our simulation. 
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Figure 3: The lepton number asymmetry, n^ = Im ( /* / j , averaged over the initial phase of 

0, as a function of time, for the same parameters as Figures ^ and |21 A lepton asymmetry 
is generated, but, as we did not include damping for the AD fields, we do not expect it to 
survive in our simulation. 



and as m^ ~ rriEw and T^ ~ 10 ^ to 10 GeV, we expect 
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we require 
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being a central value. 
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7 Conclusion 

Our baryogenesis scenario emerges from a fairly minimal extension of the MSSM with su- 
perpotential of the form 

W = XuQHuU + XdQHdd + XeLH^e + X^(l)^H^Ha + X,{LH^f + A^0^ + A^0xX- (66) 
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is a flaton field whose potential drives thermal inflation, solving the moduli problem. 
also naturally generates the /i-term of the MSSM when it acquires its intermediate scale 
vacuum expectation value, x ^^^ X represent some additional SU(5) multiplets needed to 
hold the flaton at the origin during the thermal inflation. After the thermal inflation they 
acquire intermediate scale masses from the flaton vacuum expectation value. Their effect on 
the renormalization group running of MSSM parameters from the GUT scale might provide 
a signature of our model. One could also consider embedding the Peccei-Quinn axion in the 
flaton sector. 

Our model makes the key assumption that 

"^i + <„<0 (67) 

for our AD flelds to become unstable near the end of thermal inflation. This violates a well 
known stability constraint [23 120] j implying that there exists a deep non-MSSM minimum in 
our model. This is potentially dangerous. However, we showed explicitly that the dynamics 
does not cause the flelds to become trapped in this deeper minimum. It is also known that 
the time scale for quantum tunnelling from our vacuum to this deeper minimum is larger 
than the age of our universe for a wide range of parameters [2Z] • So the above assumption 
is phenomenologically consistent and becomes a clear signature of our model which can be 
checked in future particle accelerator experiments. 

The outline of our baryogenesis scenario is as follows. We start with all flelds at the origin 
during thermal inflation. As the temperature drops, near the end of thermal inflation, the 
LHu flat direction rolls away from the origin. Then, rolls away from the origin, ending the 
thermal inflation and inducing a nonzero value for Hd- This non-zero value for Hd stabilizes 
some potentially dangerous quark and lepton fleld directions, shielding the dynamics from the 
deep non-MSSM minima discussed above. When reaches of order its vacuum expectation 
value, the back-reaction of Hd induces a lepton number violating cross-term which rotates 
the phase of LH^ generating a lepton number asymmetry. Simultaneously, L, iJ„ and Hd 
are brought back in towards the origin due to the /x-term generated by 0. We assume 
the generated asymmetry is conserved due to damping of the amplitude of the AD fleld 
oscillations by preheating, thermal friction and decay processes. While the oscillating flaton 
still dominates the energy density of the universe, the decay of the AD flelds completes, 
partially reheating the universe to a temperature high enough to restore the electroweak 
symmetry. Sphaleron processes then convert the lepton asymmetry into a baryon asymmetry. 
Finally, the flaton decay completes at a temperature of order GeV, diluting the baryon 
asymmetry to the observed value. 

As described above, despite its minimal and strongly constrained structure, our model 
has good features for baryogenesis providing rich physical content. This is a very attractive 
feature of our model. 

One key topic for future work is a proper analysis of the damping and decay of the 
AD flelds, which is crucial to the conservation of the generated lepton asymmetry. Our 
numerical simulation was only intended to model the generation of our asymmetry and 
check the stability of our ansatz, and so only simulated homogeneous fields. Thus it could 
not address the damping and decay issues. Therefore, to complete the story, a full analysis of 
the damping and decay processes and an inhomogeneous simulation including decay products 
are required. We postpone this challenging task to future work. 
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